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ABSTRACT 

We present results from X-ray and radio observations of the recently dis- 
covered young Vela-like pulsar PSR J2021+3651, which is coincident with the 
EGRET 7-ray source GeV 2020+3658. A 19.0-ks Chandra ACIS-S observa- 
tion has revealed a ~20" x 10" pulsar wind nebula that is reminiscent of the 
equatorial tori seen around some young pulsars, along with thermal emission 
from an embedded point source (/cTqo = 0.15 ± 0.02 keV). We name the nebula 
G75. 2-1-0.1. Its spectrum is well fit by an absorbed power-law model with pho- 
ton index F = 1.7lo;2, hydrogen column density uh = (7.8]'^]^;4) x lO^^cm"^, and 
an unabsorbed 0.3— 10.0 keV flux of (l.Q;!;^:^) x 10^^^ erg cm~^ s^^. We have 
spatially fit G75. 2+0.1 with a model that assumes a toroidal morphology, and 
from this we infer that the torus is highly inclined (83° ± 1°) to the line of sight. 
A 20.8-ks Chandra observation in continuous-clocking mode reveals a possible 
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pulse detection, with a pulsed fraction of ~37% and an H-test probability of 
occuring by chance of 1.2 x 10~^. Timing observations with the Arecibo radio 
telescope spanning two years show that PSR J2021+3651 glitchcd sometime be- 
tween MJDs 52616 and 52645 with parameters Az//z/ = (2.587 ± 0.002) x 10"^ 
and Az>/z> = (6.2 ± 0.3) x 10~^, similar to those of the largest glitches observed 
in the Vela pulsar. PSR J2021H-3651 is heavily scattered (r^c = 17.7 ± 0.9 ms at 
1 GHz) and exhibits a significant amount of timing noise. 

Subject headings: pulsars: general — pulsars: individual (PSR J202 1+3651) — 
stars: neutron — gamma-rays: individual (GeV 2020-1-3658) — ISM: individual 
(G75.2+0.1) 

1. Introduction 

Pulsar wind nebulae (PWNe) allow us to study the basic energetic processes of young, 
rotation-powered pulsars. It is still unclear how the rotational kinetic energy of these objects 
is converted into a relativistic particle outflow. The Chandra X-ray Observatory's unprece- 
dented spatial resolution has revealed that the winds from young, energetic pulsars can lead 
to highly structured PWNe. In some cases these PWNe take the form of an equatorial out- 
flow with polar jets, where the best example of this is the Crab nebula (Wcisskopf ct al. 
2000). These torus plus jet nebulae, of which there are roughly a half-dozen examples, are 
particularly interesting because a) they allow us to study the interaction of the PWN with its 
surrounding environment (the location of the bright torus is generally believed to mark the 
inner termination shock of the pulsar wind) and b) the well-defined morphology of the PWN 
can constrain the geometry of the pulsar. This is crucial for testing models of high-energy 
pulsations, such as the outer-gap model (e.g. Romani 1996) and the polar-cap model (e.g. 
Daugherty & Harding 1996). 

PSR J2021-I-3651 (spin period P ~103.7ms, dispersion measure DM ~369pc cm~^, 
and flux density at 1400 MHz 5*1400 ~100 /xJy) was discovered by Roberts et al. (2002) with 
the 305-m Arecibo radio telescope during a targeted search for radio pulsations from X-ray 
sources proposed as counterparts to unidentifled EGRET 7-ray sources (Roberts, Romani, & 
Kawai 2001). Timing observations of PSR J2021+3651 revealed that although it is dim, it is 
young and energetic (characteristic age = P/2P = 17kyr, E = Att'^IP/ P^ = 3.6 x 10'^^ ergs 
s~^, assuming moment of inertia / = lO''^ g cm~^), and a hkely counterpart to the ASCA X- 
ray source AX J202 1.1+3651 and the EGRET 7-ray source GeV J2020+3658. This discovery 
bolsters the idea that many of the Galactic unidentified EGRET sources (see Hartman et al. 
1999) are young energetic pulsars that were not detected in large-scale radio pulsar surveys. 
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While the spin properties of PSR J202 1+3651 are very similar to those of the known 7- 
ray pulsars Vela (PSR B0833-45, P = 89.3 ms, ^ = 6.9 x lO^^ergs s"^) and PSR B1706-44 
(P = 102.5ms, E = 3.4 x 10'^^ ergs s~^), the 12kpc distance inferred from its DM (which 
places it on the far edge of the outer spiral arm) using the NE2001 electron model of the 
Galaxy (Cordes & Lazio 2002) would make it much more 7-ray efficient than either Vela or 
PSR B1706— 44. In fact, it would have to be even more 7-ray efficient than the much-less 
energetic 7-ray pulsars Geminga and PSR B1055— 52 (assuming the NE2001 DM distance 
to PSR B1055— 52 of 0.72 kpc). This is contrary to the observation that 7-ray efficiency 
rjj = Lj/E increases with age (Thompson et al. 1999) and the theoretical expectation that 
it is inversely proportional to the square- root of the spin-down luminosity (77-y oc e.g. 
Zhang & liarding 2000). However, since DM distances can be off by a factor of a few, 
PSR J202 1-1-3651 may be significantly closer (or farther) than its DM implies. 

The similarity of PSR J202 1-1-3651 to the Vela pulsar, which is surrounded by a highly- 
structured pulsar wind nebula (see Helfand, Gotthclf, & Halpern 2001), prompted us to 
make X-ray observations of the area surrounding PSR J2021+3651. Only ^^30 of the >1300 
known pulsars in the Galaxy are Vela-like (i.e. with P ~ 100 ms, E ~ 10'^^ — 10^^ ergs s~^, 
and characteristic ages ~10— 20kyr, see Kramer et al. 2003). Half of these are known 
to have associated PWN (Kaspi, Roberts, & Harding 2004). Here we present the results of 
X-ray observations made with the Chandra ACIS-S detector. These observations resolve the 
ASCA source into a new PWN, which we name G75. 2-1-0.1, and a thermally emitting point 
source, likely the neutron-star surface. Chandra observations in continuous-clocking mode 
show a potential detection of X-ray pulsations from PSR J2021-I-3651. We also present radio 
timing observations made with the Arecibo telescope, which show that PSR J2021-I-3651 has 
glitched and is highly scattered. 



2. X-ray Observations 

2.1. X-ray Imaging 

We made Chandra observations with the ACIS-S detector in VFAINT mode on 2003 
February 12 to image PSR J2021+3651 and its surrounding region. To reduce CCD photon 
pile-up from a potentially bright point source, a quarter sub-array was used. These data 
were analyzed using CIAO^ version 3.0.1 and CALDB 2.3. Filtering for good time intervals 
and accounting for dead-time resulted in a total integration time of 19.0 ks. We subtracted 
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particle background from the image using a 100-ks combination of two ACIS-stowed images 
from 2002 Sept and 2003 May^ and then corrected for the different exposure over the chip. 
The resulting image is shown in Figure 1. 

This image reveals a point source embedded in an axisymmetric inner nebula ~20" x 10" 
across, surrounded by fainter diffuse emission. We refer to this obvious PWN as G75. 2+0.1. 
The point source is at right ascension 20''21"^05^46 and declination +36°51'04".8 (J2000). 
We estimate a positional error of ~1" by comparing the position of another bright X-ray 
point source on the S3 chip with a catalogued optical counterpart. Subtracting the nebular 
background, the point source has a count rate of 0.0094 ± 0.0009 counts s~^ and the nebula 
has a background-subtracted count rate of 0.075±0. 004 counts s~^ in the 0.3— 10.0 keV range. 

The greyscale of Figure 1 shows the image smoothed with a 1.2" (FWHM) gaussian to 
bring out the fine structure of the inner nebula. The white contours indicate the position of 
the point source and are shown to distinguish it from the bright bar running along the axis 
of the nebula. This bar is the most obvious resolved feature of the nebula. Also overlaid 
are black contours of the same image smoothed with a 3.5" gaussian and scaled to bring out 
what are possibly faint jets along the axis of the nebula. We estimate 62±32 counts may be 
coming from these putative jets. The nebula appears symmetric about the minor axis, but 
the point source is clearly offset from the major axis of the nebula. The overall morphology 
is suggestive of the torus plus weak jet morphologies of the Crab and (likely) Vela nebulae 
(Weisskopf et al. 2000; Helfand et al. 2001). We discuss the geometry of the PWN further 
in §4.4. 



2.2. X-ray Spectroscopy 

The spectra of the PWN and point source can be measured individually using Chandra 's 
excellent spatial resolution. We fit the spectrum of the nebula and point source using ClAO's 
jS/ierpa environment (version 3.0.1) and XSPEC (version 11.2.0). To isolate the nebula, we 
used a circular annulus centered on the point source and extending from a radius of 1.6—100 
pixels (0.492" x 0.492" per pixel) and a circular region from an apparently fiux-free area 
of the S3 chip to estimate the background. We find that the nebula has a hard spectrum 
(see Figure 2) that is well fit (reduced Xi^'^ — 1-1 for 27 degrees of freedom) by an absorbed 
power-law model with hydrogen column density uh = (7-8lJ;4) x 10^^ cm~^, photon index 
r = 1.7toi, and a 0.3-10 keV absorbed fiux of (1.2 ± 0.1) x lO'^^ ^^-2 g-i 
confidence intervals). The best-fit black-body model gives reduced x,^^ = 2.3 for 27 degrees 
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of freedom, a comparatively poor fit to the data. 

To isolate the point source, we used a circular extraction region of radius 1.6 pixels, and 
a surrounding annulus of width 1.5 pixels to subtract the nebular background. Because the 
energy dependent PSF extends significantly past 1.6 pixels at this location, we corrected for 
this by using "mkpsf to generate model PSFs at several energies between 0.5— 6keV from 
which we estimated the source minus background enclosed energy fraction. We then fit a 
linear function to these points which we used to correct the effective area model used in 
the spectral fit. Based on the observed count rate, we estimate a negligible pile- up fraction 
of ~0.6% using HEASARC's wcbPIMMs tool^. Wc find that there is very little emission 
from the point source above 3keV. Its spectrum is adequately fit (reduced Xu^ = 1-7 for 
9 degrees of freedom) by an absorbed black-body model, where we have fixed nn to the 
value obtained by fitting the PWN's spectrum. The best-fit black-body model and 90% 
confidence region is kT^o = 0.15 ± 0.02 keV with an absorbed fiux of (2.8 ± 0.2) x 10~^^ erg 
cm~^ s~^ (0.3— 10.0 keV). Since kT^o and uh arc highly covariant, we also made a joint fit of 
an absorbed black body for the pulsar and an absorbed power law for the nebula, tying the 
uh of both models together. This gave values and errors that were essentially the same as 
those from the individual fits. Since wc expect non-thermal emission from the magnetosphere 
of the pulsar, we also fit a black-body plus power-law model to the point source, fixing the 
photon index at a typical value of 1.5. This did not improved the fit (reduced x^"^ = 1.8 for 
7 degrees of freedom) , but the temperature and fiux of the black-body component were left 
relatively unchanged. This suggests a maximum non-thermal absorbed fiux of ~ 3 x 10~^^ 
erg cm-2 (0.3-10.0 keV). 

The total flux measured here is a factor of four smaller than seen by ASCA (Roberts 
et al. 2001; Roberts et al. 2002). We note that the ASCA extraction region has a radius 
of ~4', and thus includes areas well outside the 2'- wide field of view of the quarter array 
used here. A highly smoothed image of the Chandra data suggests there may be additional 
diffuse emission that could extend well beyond the chip. Furthermore, examination of the 
ASCA image suggests there may be additional point sources or structures within the ASCA 
extraction region that are not within the Chandra field of view. However, despite the fiux 
discrepancy, the derived uh, nebular photon index and point source kT^ are consistent with 
the ASCA measurements of the source. 
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2.3. X-ray Pulsations 

To search for X-ray pulsations from PSR J2021+3651, we also obtained Chandra ACIS- 
S data in continuous-clocking mode, which provides an effective time resolution of 2.85 ms 
by sacrificing one dimension of spatial resolution. The roll angle of the telescope placed 
the 1 dimensional image at an angle of ~62° (N-E). The events were filtered for good time 
intervals and periods of high background, leaving 20.8 ks of integration time. We corrected 
the read-out times to intrumental times of arrival by removing telescope dither, motion of 
the detector with respect to the telescope, and the average lag between the readout time 
and the true times of arrival of the photons. The events were then barycentered using the 
JPL DE405 ephemeris (Standish 1998). Given prior knowledge of the point-source spectrum 
from our imaging observations, we produced an event list with a 0.5— 3keV energy cut and a 
3-pixel-radius extraction region. We also produced event lists for a number of other energy 
cuts and a range of extraction regions and then dithered the events by adding random 
deviates in the range 0— 2.85 ms to remove windowing effects. The event hsts were folded 
at the predicted spin frequency of 9.64094861 Hz according to the contemporaneous radio 
timing ephemeris (see §3), which revealed a hkely periodic signal (see Figure 3). The signal 
shows up for a range of extraction region sizes and energy cuts, but has a maximum signal 
strength for an extraction region of radius 3 pixels with a 0.5— 3keV energy cut. Performing 
a variant of the H-test (De Jager et al. 1989), using unbinned events, we find an H-test score 
of h ~11.4, which corresponds to a chance probability of 1.2 x 10~^ for a single trial. This 
has an equivalent Gaussian sigma of 3.7a. Searching to a maximum harmonic of 20, the 
H-test finds the largest significance when two harmonics are summed. 

Although this signal is not a certain detection of pulsations from PSR J2021+3651, it 
is intriguing. If real, these pulsations indicate a pulsed fraction of ~37%, where we define 
pulse fraction to be the ratio of the number of counts above the 1-a upper limit on the lowest 
bin in the fight curve (see Figure 3) to the total number of background-subtracted counts 
(in this case, most of the background is from the nebula). From the point source count-rate 
of our imaging observations, we estimate that ~ 190 of the counts in the light curve are from 
the pulsar and that the rest are from the nebula and other background. We discuss the pulse 
shape in §4.2. 

3. Radio Timing 

Since the discovery of PSR J202 1-1-3651 in 2002 February, we have been making roughly 
monthly radio timing observations of the pulsar at 20 cm using the Arecibo radio telescope 
and multiple Wideband Arecibo Pulsar Processor digital correlator backends (WAPPs, for 
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details see Dowd, Sisk, & Hagen 2000). Each WAPP was configured for 100 MHz of band- 
width, with 512 lags and 200 fis samples. Whenever possible, we used three WAPP backends, 
typically centered at 1170, 1420, and 1520 MHz. PSR J2021+3651 is significantly scattered 
(see Figure 4). Fitting the profiles with a gaussian convolved with a one-sided exponen- 
tial (e"*/^*'') we find a frequency dependence of Tgc that is well fit by a power law of index 
—4.2 ±0.6. This is consistent with a pure Komolgorov spectrum (i.e. r^c oc u^'^-'^, see Rickett 
1977). Using the measured frequency dependence, we find Tgc — 17.7 ± 0.9 ms at 1 GHz (90% 
confidence) . 

Pulse profiles from each epoch were cross-correlated with a representative high signal- 
to- noise template to generate times of arrival (TOAs). Using TEMPO^ we fit these TOAs 
to a multi-parameter rotational model of the pulsar. The position of the pulsar was fixed 
at the Chandra X-ray position. Our timing analysis reveals that a glitch occured between 
MJDs 52616 and 52645. We can connect phase across it by fitting an instantaneous jump 
in frequency and frequency derivative (Av/i' — 2.587 ± 0.002 x 10~^ and Az>/z> = 6.2 ± 
0.3 X 10^'^), which are comparable to those seen in the largest glitches of the Vela pulsar 
(Dodson et al. 2002). There is only marginal evidence for a glitch decay in the data, which 
is not surprising given that the typical timescale for such a decay is on the order of days to 
weeks (Dodson et al. 2002) and is therefore not resolved by our sparse monitoring. The pre- 
and post-ghtch rotational ephemerides of PSR J202 1-1-3651 and the complete set of glitch 
parameters are presented in Table 1, along with derived quantities. 

The post-fit timing residuals show significant systematics when only the frequency and 
first frequency derivative are fit (see Figure 5). We find that the residuals can be rendered 
featureless by including a frequency second and third derivative into the fit. Alternately, the 
residuals can be rendered featureless by allowing TEMPO to fit for position. However, the pre-, 
post-glitch, and Chandra positions are inconsistent, suggesting that timing noise dominates 
the residuals rather than an incorrect position. 

4. Discussion 
4.1. Distance 

As discussed in Roberts et al. (2002), the DM distance of PSR J202H-3651 is quite 
large (~12 kpc using the NE2001 electron density model of the galaxy, see Cordes & Lazio 
2002) and would place PSR J2021-I-3651 on the outer edge of the outer spiral arm of the 
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Galaxy. This large distance leads to a very high inferred 7-ray efficiency {rjj = L^/E ~ 0.18 
in the 100 MeV— 10 GeV range, assuming 1-sr beaming). Cordes & Lazio (2002) note that 
the electron density model is seldom off by more than 50% when predicting distances, except 
when the line of sight intersects a particularly dense scattering screen. Wc find a scattering 
measure SM = (1.5 ± 0.1) d^^^^pc m-^o/s towards PSR J2021+3651, where Tsc = 1.10ms 
SM^/^i/~^^/^(i (Cordes & Lazio 2002), that is two orders of magnitude larger than predicted 
by the NE2001 model (although PSR J2021+3651 is in the direction of the Cygnus region 
and pulsars in that direction are often highly scattered). This suggests the presence of a 
clump of scattering material along the line of sight not accounted for by the NE2001 model. 
However, small clumps of material that only change the DM by 10— 20pc cm~^ (< 5% in 
the case of PSR .12021+3651) can change the SM by a factor of a hundred. Hence, the 
large SM does not by itself imply a significantly closer distance. In order for the distance 
to imply a much smaller 7-ray efficiency, a significant source of dispersion (~200pc cm~^) 
unaccounted for by the model is necessary. Such a DM enhancement could come from an 
ultra-compact Hll region or an OB association. An examination of the Midcourse Space 
Experiment 8.3 //m IR data (available from the NASA/IPAC Infrared Science Archive) and 
preliminary radio imaging with VLA data we have obtained (work in preparation) show no 
evidence for an ultra-compact Hll region. The OB association Cyg 1,8,9 lies in the direction 
of PSR J2021-F3651(Yadigaroglu & Romani 1997; Mel'nik & Efremov 1995), however there is 
currently no evidence that this is significantly contributing to the amount of ionized gas along 
the line of sight or that the pulsar was born in this association. Since the DM distance places 
PSR J2021+3651 at the very outer edge of the Galaxy one would expect our spectral fit to 
give an uh consistent with the total column density in this direction. The X-ray absorption 
found here, ur = (T.S^J;!) x lO^^ cm ^, is lower than the total Galactic Hi column density in 
the direction of PSR J2021+3651, 1.2 x 10^^ cm-^ estimated using the FTOOL nh^ (which 
uses the Hi map of Dickey & Lockman 1990). This argues for a distance somewhat closer 
than the predicted DM distance, but not closer by factor of a few. Hence a distance of 
~8kpc is quite plausible, but a distance of less than 5kpc, as required to make the 7-ray 
efficiency similar to that of PSR B1706— 44 (j]^ = 0.01, sec Roberts 2003), is hard to justify 
given current observations of the region. We adopt a nominal distance of lOkpc to use in 
scaling quantities {dio = d/lOkpc) that depend on distance. 
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4.2. Neutron Stcir Thermal Emission and Pulsations 

We find a black-body temperature of kT^o = 0.15 ± 0.02 keV for the point source, 
similar to that found for the Vela pulsar {kToo = 0.128 ± 0.003 keV; Pavlov et al. 2001) and 
PSR B1706-44 {kT^ = 0.14 ± 0.02 keV; Gotthelf et al. 2002). Accounting for the effects 
of gravitational redshift, this corresponds to a surface temperature kTg — 0.19 ± 0.03 keV 
for the canonical neutron-star radius of 10 km and mass of 1.4 M©. The point source has 
an unabsorbed flux Fx = 6.3 x 10~^^ erg cm~^ s~^ in the range 0.3— 10.0 keV. Scahng the 
luminosity to a distance of lOkpc, Lx = 7.5 dio^ x 10'^'^ erg s~^ (0.3— 10.0 keV). One can 
calculate the ratio of the black-body emission radius -Re and the distance to the source using 
the bolometric flux 7.8 x 10^^'^ erg cm~^ s~^ and the Stefan-Boltzman law. Correcting for 
gravitational redshift, this gives Rio/dio = 0.3, where Riq = i?e/10km. For any reasonable 
distance, the simple black-body model suggests that the emission is coming from hot spots 
on the surface, rather than from the entire surface. This is not unreasonable, as some 
high-energy emission models predict hot spots due to back-flowing particles. 

Using the non-magnetized neutron-star hydrogen and iron atmosphere models of Gaen- 
sicke et al. (2002), which assume a neutron-star radius of 10 km and a mass of 1.4 Mq, changes 
the temperature signiflcantly. Fixing the uh to the value derived from the nebular flt, for a 
hydrogen atmosphere we get a best-flt value of kTg — 0.06 ± 0.01 keV (which is similar to the 
hydrogen atmosphere flt of the Vela pulsar, kT^ — 0.059 ± 0.003 keV; Pavlov et al. 2001) 
and a flux of 2.9 x 10~^^ erg cm~^ s~^ (0.3— 10.0 keV). However, the normalization from this 
flt implies a ~0.8kpc distance, which is unreasonable. Again flxing the uh to the nebular 
value, for an iron atmosphere model we get a best-fit value of kTs = 0.22 ± 0.01 keV and a 
fiux of 3.6 x 10^^^ erg cm~^ s~^ (0.3— 10.0 keV). The implied distance here is ~13kpc, which 
is roughly consistent with the emission coming from the entire surface. 

Based on the soft spectrum of the point source, one would expect pulsations to be 
due primarily to thermal emission. The profile found by folding the X-ray events of our 
continuous-clocking data with the radio ephemeris shows a somewhat narrow (duty cycle 
~0.3) peak and interpulse, which is possible for thermal pulsations given the right geometry. 
If we are seeing emission from hot spots at both poles then it is possible to get two sinusoidal 
peaks separated by a plateau of emission (Beloborodov 2002). In fact, the geometry inferred 
from the PWN (see §4.4) is consistent with the pulse profile seen here. As noted in §2.2, 
a non-thermal component is also allowed by the spectrum. The possible pulsations might 
be a combination of pulsed thermal and magnetospheric emission, but the number of pulsed 
counts (~69) is too large to consist solely of magnetospheric emission. All the known 7-ray 
pulsars have also been detected as X-ray pulsars with non-thermal emission components to 
their spectra (Kaspi et al. 2004). We note that a lack of X-ray pulsations would not preclude 
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PSR J2021+3651 from being a 7-ray pulsar. For instance, X-ray pulsations in PSR B1706— 44 
are relatively difficult to detect even though PSR B1706— 44 is a well-established 7-ray pulsar. 
Thus, PSR J202 1-1-3651 remains an excellent 7-ray pulsar candidate. 

4.3. PWN Spectrum 

We found that the spectrum of the nebula is well fit by an absorbed power-law model 
with uh = i^-Stl'l) X 10^^ cm~^ and F = l-T^oi- The photon index F is roughly consistent 
with what is observed from the PWN of other similar pulsars like Vela (F = 1.50±0.04; ) and 
PSR B1706-44 (F = 1.34+°J^; Gotthelf et al. 2002), although it is somewhat higher. The 
relationship of Gotthelf (2003) between the photon index of the nebula and the spin-down 
energy of the pulsar (equation 3 of that paper) gives F = 1.22 ± 0.33, which is marginally 
consistent with the F measured here. 

The unabsorbed nebular flux of G75. 2+0.1 is Fx = (LQlgg) x 10^^^ erg cm^^ s^^ in the 
0.3— 10.0 keV range. This corresponds to an X-ray luminosity in the same energy range of 
Lx = (2.3 dio^) X 10^^ erg s^^. For the same energy range, we calculate an X-ray efficiency 
(r/x = Lx/E) of 0.009 dio^ from the combined point-source and nebular luminosities. This 
value fits well with the claimed relationship between the spin-down energy of a pulsar and 
its nebular X-ray efficiency (c.f. Possenti et al. 2002; Chevaher 2000) even if the distance is 
as close as 3kpc. However, note that this is only for the bright inner nebula. If there is 
significantly more flux coming from diffuse emission (as suggested by the ASCA data and in 
analogy with the Vela nebula, see below) then the efficiency will be higher. 

4.4. PWN Size and Morphology 

The generally accepted view of young PWNe (see Kaspi et al. 2004, and references 
therein) is that they are synchrotron bubbles blown at the center of an expanding super- 
nova remnant (SNR). High-energy electrons and positrons, along with magnetic field, are 
continuously injected into the bubble by the pulsar. Given that the characteristic age of 
PSR J2021+3651 is much larger than the typical time scale for a SNR entering the Sedov 
phase, it is likely (although not certain since characteristic ages can be off by an order of 
magnitude or more) that the outer edge of the bubble (at radius Rp) is expanding subsoni- 
cally into the shocked supernova ejecta. The bubble is bounded at the inner edge (at radius 
Rt) by the termination of the ultra-relativistic pulsar wind. 

The physical size of G75.2+0.1's inner nebula is - (0.5 x 1) dw pc. As PSR J2021+3651's 
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spin-down energy is within a factor of two of that of the Vela pulsar, one might naively expect 
the sizes and morphologies of their nebulae to be comparable. Scaling the ~ 10" radius of 
G75.2+0.1's inner nebula to that of the Vela pulsar's inner nebula, which has a radius of 
52. "7 (Helfand et al. 2001) and whose distance is well known (290 pc, Dodson et al. 2003), 
places G75. 2+0.1 at a distance of ~1.5kpc, which seems unreasonable given the DM. In fact, 
G75. 2+0.1 would have to be a few times larger than the Vela nebula to put it at a distance 
that seems at all reasonable based on its DM. However, since the Vela inner nebula is only 
the central part of a much larger and more luminous nebula, a better comparision between 
the sizes of these two objects could be made if the true extent of G75. 2+0.1 is determined 
in a deeper observation with a larger field of view. 

The bright inner nebula of Vela and G75. 2+0.1 may be directly downstream of the 
termination shock of the relativistic pulsar wind. The location of the termination shock 
depends primarily on the energy output of the pulsar and the pressure in the nebula. The 
pressure in the nebula depends upon the total energy input during the pulsar's lifetime and 
the losses due to adiabatic and synchrotron processes (possibly complicated by the passage 
of the SNR reverse shock). The total energy input is highly dependent on the initial spin 
period, which is unknown and could range from ~10— 100 ms (Kaspi et al. 2001; Migliazzo 
et al. 2002). We currently do not know the full extent of G75. 2+0.1, which is likely much 
larger than the inner nebula seen here. While it is true that the pressure in the PWN 
should be balanced by the pressure in the surrounding SNR, we know neither the age (as the 
characteristic age is a poor estimate of the true age for young pulsars) nor the size (since we 
have not observed a shell) of the potential remnant. Therefore, although the Vela pulsar is 
twice as energetic as PSR J202 1+3651, it is entirely possible that G75. 2+0.1 is several times 
larger than the Vela nebula. 

The size of the torus places an upper limit on the termination shock radius Rt- In 
analogy to Eqn. 1 in Roberts et al. (2003) (see also Chevaher 2000) we can estimate the 
post-shock magnetic field to be 

Sjv~2 X 10-H ft'''i\^ )G 

^ ft(RT /r)dio ^ 

where r is the measured location of the torus, < 1 is the ratio of the magnetic energy to 
the total energy in the post-shock flow and /t < 1 is a geometrical factor allowing for non- 
spherical outflow. Note that the PWN of PSR B1706— 44 (which has a virtually identical 
E to PSR J2021+3651) is both much fainter and smaller than G75.2+0.1. This implies 
that the magnetic field in PSR B1706— 44's nebula is smaller but the density of synchrotron 
emitting particles is larger. PSR B1706— 44 is a known TeV source (Kifune et al. 1995), this 
suggests PSR J202 1+3651 may be an excellent target for ground-based Cerenkov detectors. 
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Our G75. 2+0.1 image is rather sparse and Poisson- limited, and so a unique model for its 
morphology cannot be defined. Nevertheless, the axisymmetry of G75. 2+0.1 suggests that 
the morphology may be best described as that of a torus, with possibly weak jets directed 
along the axis of symmetry (see Figure 1). To investigate the question of G75.2+0.1's mor- 
phology quantitatively, we spatially fit the PWN using the relativistic torus fitting described 
in Ng & Romani (2004). Since G75. 2+0.1 appears to have a double ridge, our fiducial model 
is similar to that of the Vela pulsar PWN: twin tori of radius r, symmetrically spaced by 
distance d on either side of the pulsar spin equator. The torus axis is at position angle ^ 
(N-E) and inchnation ( to the line of sight. The bulk flow velocity down-stream from the 
termination shock is (3t and a 'blur' parameter 5 is applied to account for the unmodeled 
post-shock synchrotron cooling length. A uniform background and a detailed point source 
model computed for the pulsar spectrum and chip position are also included. The 3-D model 
is projected to the plane of the sky, a Poisson-based figure of merit is minimized to find the 
best-fit parameters and the statistical errors on these parameters are estimated by refitting 
many Poisson realizations of the best-fit model, each with the same number of counts as 
seen in the data. 

The derived parameters from the fit are summarized in Table 2. The errors quoted 
are statistical l-cr values only. Systematic errors due to unmodeled features are difficult 
to quantify. One estimate of their magnitude can be made by excluding an elliptical zone 
about the pulsar, extended along the symmetry axis where an unmodeled polar jet might 
be present. Re-fitting the model to the surrounding data, we obtain values within 1.5 a 
of the original fit for all parameters, except Pt which now is 0.64 ± 0.02; Ng & Romani 
(2004) note that for this parameter the statistical errors tend to be small and dominated 
by systematic uncertainties for edge-on tori. In this particular case, the value of /3t is likely 
being dominated by the bright clump to the south-west of the pulsar. We have also fit 
a single torus; again the geometrical parameters are quite similar, although a larger blur 
parameter 6 = 1.8" is required to cover the broader torus. The fit statistic for the model is 
of course larger than that of the 2-torus model, by an amount similar to that corresponding 
to l-(7 excursions in the individual fit parameters. Unfortunately, this does not alone show 
that the double torus is a significantly better model; we cannot apply an F-test, since the 
fit statistic is not x^. 

Clearly a deeper image is needed to draw firm conclusions, but we can discuss the 
geometrical implications of the fit angles. The large C implies that we are viewing the pulsar 
near the spin equator. In the outer magnetosphere picture of Romani Sz Yadigaroglu (1995) 
this is indeed where we expect to see 7-ray emission. More particularly. Figure 4 of that 
paper shows that for ( f» 83° and magnetic impact angle j3 modest (to allow detection of the 
radio emission), we expect two 7-ray pulses separated by ~180°, with the first lagging the 
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radio peak by ~10°. Intriguingly, McLaughlin & Cordes (2003) found apparently significant 
pulsations folding the 7-ray photons from this source near the extrapolated P for two viewing 
periods; the light curves show two narrow peaks separated by 0.5 in phase. It is important 
to note that for ( near 90°, such a light curve is also naturally expected in the two pole 
model of Dyks & Rudak (2003) and may be possible in polar-cap models (Daugherty & 
Harding 1996). The separation of the two tori suggest bright emission at spin axis angles 
9t — tan~-^(2rsinC/(i) = 77°, 103°, as viewed from the pulsar. If these represent enhanced 
e^ flowing radially from near the magnetic poles then we might infer a magnetic inclination 
a ~ 77°. This would imply an impact parameter (3 — ( — a ^ +6°. Polarization sweep 
measurements could test this picture. The unscattered pulse half-power width ~33° of 
PSR J2021-I-3651 is substantially wider than the 2.5°(P/1 s)~^/^/sinQ; width expected for a 
core component, but is similar to the 5.8°(P/1 s)^^/^ width for an outer cone (Rankin 1993). 
Of all the fit geometrical parameters, the most robust is the fit PA of the symmetry axis 
Although difficult to measure, given the large DM and expected substantial Faraday 
rotation, the absolute position angle at the polarization sweep maximum would provide an 
interesting check of the inference that this represents the pulsar spin axis. 



5. Conclusions and Further Work 

We have used the Chandra X-ray Observatory to observe the young and energetic pulsar 
PSR J2021-I-3651 and have detected a new PWN, G75. 2-1-0.1, as well as an embedded point 
source near the radio timing position of the pulsar. The morphology of G 75. 2-1-0.1 is perhaps 
best decribed as an equatorial torus, and we have derived geometrical parameters for the 
orientation of the system by spatially fitting the nebula. Spectral fitting of the nebula shows 
that it is similar to the nebulae of other energetically-similar pulsars. Spectral fitting of 
the point source shows that it is well fit by an absorbed black-body model or a black-body 
plus power-law model. Radio timing observations of PSR J2021-I-3651 have revealed a glitch 
similar to the largest glitches seen in the Vela pulsar. Analysis of count inuous- clocking data 
from Chandra indicates that there may be weak X-ray pulsations from the source (pulsed 
fraction of ~37%) at the period predicted by the radio ephemeris. PSR J2021-I-3651 remains 
the most likely counterpart to the high-energy EGRET 7-ray source GeV J2020-I-3651 (which 
is the 10th brightest 7-ray source above 1 GeV). It is an excellent 7-ray pulsar candidate for 
the AGILE and GLAST missions, although a contemporaneous ephemeris will be essential, 
given the large amount of timing noise the pulsar exhibits. 

Deeper observations with Chandra and XMM-Newton are needed to clarify the morphol- 
ogy of G75. 2-1-0.1 and to confirm the presence of pulsations. Such observations are also able 



-14- 



to confirm whether the faint outer jet hinted at here is real. Currently, there only about a 
half dozen PWNe with torus plus jet morphologies, with the Crab and (likely) Vela nebulae 
as the best examples. With a deeper observation, G75. 2+0.1 may prove to be the third 
best example. Further work by our group will also include analysis of radio polarization 
data we have taken with Arecibo. These data may help bolster the geometrical interpre- 
tation described here if the swing of the polarization angle across the pulse is measurable. 
Proper- motion measurements could check if PSR J202 1-1-3651 is moving along the axis of 
the nebula, which is the case for both Vela and the Crab. However, measuring the proper 
motion of PSR J202 1-1-3651 (which will hkely be small since the distance is hkely large) will 
be very difficult both through radio timing because of timing noise or through VLBI because 
the radio source is very faint. We are also analyzing VLA observations (in A, C, and D 
arrays) of the source region to look for a radio PWN and/or SNR. This is important for 
characterizing the surrounding medium, which has an important effect on the morphology 
of the PWN and its size. 
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Table 1. Measured and Derived Parameters for PSR J202 1+3651 



Parameter 



Value 



Right ascension (J2000.0)" 
Declination (J2000.0)^ 
Galactic longitude / (deg)' 
Galactic latitude b (deg)^ 
DM (pc cm~^) 

Pulse frequency u (s~^) 
Frequency derivative z> 
Epoch (MJD) 
Range of MJDs 
RMS Residuals 



Preglitch 



ms 



20''21'"05^46(5) 

+36°51'04".8(7) 

75.23 

+0.11 

369.2(2) 

9.64113518458(16) 

-8.886500(24) x lO'^^ 

52407.389 

52305-52616 

3.0 



Pulse frequency (s ^) 
Frequency derivative z> (s~^) 
Epoch (MJD) 
Range of MJDs 
RMS Residuals 



ms 



Epoch (MJD) 



Postglitch 

9.64091197572(12) 
-8.932592(12) x 10"^^ 
52730 

52645-53009 
9.3 

Glitch Parameters 

52629.97-52630.16 



Frequency jump Au/u 2.587(2) x 10~^ 

Frequency derivative jump Au/i/ 6.2(3) x 10~^ 

Derived Parameters 
Spin-down Luminosity (ergs s~^) 3.4 x 10^^ 
Surface Magnetic Field (G) 3.2 x 10^^ 

Characteristic Age = P/2P (kyr) 17 



Note. — Figures in parentheses represent uncertainty in the 
least-significant digits quoted equal to 3 times errors given by 
TEMPO. 
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'^Coordinates of the Chandra X-ray point source. We estimate 
a positional error of ~1" by comparing the position of another 
bright X-ray point source on the S3 chip with a catalogued optical 
counterpart. 

""E = An'^IP/P^ with I = 10^^ g cm^. 

'^Assuming standard magnetic dipole spindown: B = 3.2 x 
10i9(pp)i/2 Gauss. 



-19- 



Table 2. Torus Fit Parameters 



Parameter^ 


Value 




47° ± 1° 


c 


83° ± 1° 


r 


8.6" ± 0.3" 


5^ 


1.2" 




~0.8 


d 


3.7" ± 0.2" 



^Explanation of parame- 
ters in §4.4 

'^Held fixed in the global 
fit. 

*^The Doppler boosting 
is likely being dominated 
by the bright clump to the 
south-west of the pulsar. 
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Fig. 1.— A 19.0-ks Chandra ACIS-S image (0.3-7.0 keV) of PSR J2021+3651 and its inner 
nebula. The greyscale shows the image smoothed with a 1.2" (FWHM) gaussian and scaled 
to bring out fine structure in the nebula. The black contours are the same image overlaid, 
but smoothed with a 3.5" gaussian and scaled to bring out faint extended emission including 
what may possibly be faint jets along the axis of the nebula. The white contours indicate 
the point source and are shown to distinguish it from the bright bar running along the axis 
of the nebula. The dashed fine shows the best-fit position angle ^ = 47° ± 1° of the torus 
fit (sec §4.4), which in analogy with the Crab pulsar and nebula likely corresponds to the 
projected spin axis of the pulsar. 
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Fig. 2. — Left: The background-subtracted spectrum of the point source, fit using an ab- 
sorbed black-body model in Sherpa (CIAO 3.01). The best-fit temperature and 90% confi- 
dence interval is kT^ = 0.15±0.02keV with a 0.3-10.0 keV absorbed flux of (2.8±0.2) xlO"^^ 
erg cm~^ s~^. The uh was fixed at 7.8 x 10^^ cm~^, that found by fitting the nebular spectrum. 
Right: The background-subtracted spectrum of the nebula fit with an absorbed power-law 
model. The best-fit parameters and 90% confidence intervals are uh = (T-^-ia) ^ 10^^ cm~^, 
r = 1.7toL and an absorbed flux of (1.2 ± 0.1) x 10"^^ erg cm-^ s"^ (0.3-10.0 keV). 
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Fig. 3. — Chandra continuous clocking data of PSR J2021+3651 folded using a local radio 

epliemeris for the pulsar. Two cycles are plotted for clarity and there are ten bins across the 
profile. The phase of the radio pulse is indicated by an arrow. The dashed line indicates the 
approximate background level. Based on this, there are ~69 pulsed counts and we estimate 
that ~190 counts are from the pulsar. This gives a pulsed fraction of ~37%. 
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Fig. 4. — Summed pulse profiles of PSR J2021+3651 at center frequencies of 1175 MHz (top) 
and 1475 MHz (bottom). The peak fluxes of the two proflles are normalized and the pulse 
phases aligned at the pulse peak. By fitting the profiles to a gaussian convolved with an 
exponential we calculate a scattering time Tsc — 17.7 ± 0.9 ms at 1 GHz. 
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Fig. 5. — Residuals from phase-coherent timing of PSR J2021+3651 fitting only for period 
and the first period derivative. The residuals can be rendered featureless by fitting for a 
period second derivative. Alternately, the residuals can be rendered featureless by fitting for 
position (instead of keeping it fixed at the X-ray position). However, the fitted positions from 
the pre- and post-fit ephemerides are inconsistent suggesting that timing noise dominates 
the residuals, rather than an incorrect position. Left: pre-glitch residuals. Right: post-glitch 
residuals. 




Fig. 6.— Left: 32" x 32" cutout of the central PWN. Middle: best-fit double torus plus point 
source model. Right: best-fit single torus plus point source model. 



